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ABSTRACT 



We describe a method to measure the magnetic field orientation of coronal 
mass ejections (CMEs) using Faraday rotation (FR). Two basic FR profiles, 
Gaussian-shaped with a single polarity or "N"-like with polarity reversals, are 
produced by a radio source occulted by a moving flux rope depending on its 
orientation. These curves are consistent with the Helios observations, providing 
evidence for the flux-rope geometry of CMEs. Many background radio sources 
can map CMEs in FR onto the sky. We demonstrate with a simple flux rope that 
the magnetic field orientation and helicity of the flux rope can be determined 2-3 
days before it reaches Earth, which is of crucial importance for space weather 
forecasting. An FR calculation based on global magnetohydro dynamic (MHD) 
simulations of CMEs in a background heliosphere shows that FR mapping can 
also resolve a CME geometry curved back to the Sun. We discuss implementation 
of the method using data from the Mileura Widefield Array (MWA) . 

Subject headings: Faraday rotation — magnetic fields - - Sun: coronal mass 
ejections 



1. Introduction 



Coronal mass ejections (CMEs) are recognized as primary drivers of interplanetary 
disturbances. The ejected materials are often associated with large southward magnetic 
fields which can reconnect with geomagnetic fi elds and produce storms in the terrestrial 
environment (jDungeylll96ll ; iGosling et al.lll99ll ). Determination of the CME magnetic field 
orientation is thus of crucial importance for space weather forecasting. However, nearly all 
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atoms are ionized at the coronal temperature ~ 2 x 10 6 K, making it difficult to detect 
the coronal magnetic field through Zeeman splitting of spectral lines as is routinely done 
for the photospheric field. A typical way to estimate the coronal magnetic field above 2 
Rq (Rq being the solar radius) is theoretical extrapolation using the photospheric fields 
as boundary conditions, which can only be checked by comparison to the field strength 
measured from radio bursts and the orientation determined from soft X-ray observations. 
The field orientation is also hard to infer from white-light coronagraph images. Spacecraft 
near the first Lagrangian point (LI) mea sure the local fie l ds but can only giy e a warning 
time for arrival at Earth of ~ 30 minutes (jVogt et al.ll2006l ; IWeimer et al.ll2002l ). 



A possible method to measure the coronal magnetic field is Faraday rotation (FR), the 
rotation of the polarization plane of a radio wave as it traverses a magnetized plasma. The 
first FR experime nt was conducted in 1968 by Pioneer 6 during its superior solar conjunction 
( jLevy et al.lll969l ). The observed FR curve features a "W" -shaped profile over a time period 
of 2-3 hours with rotation angles up to 40° from the quiesc ent baselin e. Th is FR event was 
i nterp reted as a coronal streamer stalk of angular size 1-2° (jWoolll997l ). but iPatzold fc Bird 
( 119981 ) argue that the FR curve is produced by the passage of a series of CMEs. Joint coro- 
nagraph observations are needed to determine whether an FR transient is caused by CMEs. 
Subsequent FR observations b y the Pioneer and Helios spacecraft reveal important informa- 



tion on t he quiet coronal field (IStelzried et al 



tuations (IHoUweg et al.l Il982l ; lEfimov et al. 
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1996 



Patzold et al.lll987) and magnetic fluc- 



Andreev et al. 1997: Chashei et al. 1999 



20001 ) . FR fluctuations are currently the only source of information for the coronal field 
fluctuations. Independent knowledge of the electron density, however, is needed in order to 
study the background field and fluctuations. 

Joint coronagraph and FR measurements of CMEs were also conducted when the Helios 
spacecraft, with a downlink signal at a wavelength A = 13 cm, was occulted by CME plasma. 
Bird et al.l (119851 ) establish a one-to-one correspondence between the SOLWIND white-light 
transients and FR disturbances for 5 CMEs. Figure 1 displays the time histories of FR 
and spectral broadening for two CMEs. Note that the spectral broadening is proportional 
to the plasma density fluctuations; the increased spectral broadenin g is consistent wi th the 
enhanced density fluctuations within CMEs and their sheath regions (ILiu et al.ll2006bl ). The 
FR through the 23 October 1979 CME shows a curve (note a data gap) which seems not to 
change sign during the CME passage; a single sign in FR indicates a monopolar magnetic 
field. The 24 October 1979 CME displays an FR curve which is roughly "N"-like across 



the ze ro rotation angle, indicative of a dipolar field. Other CMEs in the work of iBird et al. 
( 119851 ) give similar FR curves, either an "N"-type or a waved shape around the zero level. 



Based on a simple slab model for CMEs, the mean transient field magnitude is estimated to 
be 10 - 100 niG scaled to 2.5i? Q , which seems larger than the mean background field. The 
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CME field geometry, as implied by these FR curves, will be discussed below. These features 
demonstrate why radio occultation measurements are effective in detecting CMEs. 

FR experiments using natural radio sources, such as pulsars a nd quasars, have a lso been 
performed. FR observations of this class were first conducted by iBird et al.l (119801 ) during 
the solar occultation of a pulsar. The advantage of using natural radio sources is that many 
of these sources are present in the vicinity of the Sun and provide multiple lines of sight 
which can be simultaneously probed by a radio array. We can thus make a two-dimensional 
(2-D) mapping of the solar corona and the inner heliosphere with an extended distribution 
of background radio sources. 

In this paper, we show a method to determine the magnetic field orientation of CMEs 
using FR. This method enables us to acquire the field orientation 2-3 days before CMEs 
reach Earth, which will greatly improve our ability to forecast space weather. The data 
needed t o implement this technique will be available from the Mileura Widefield Array 
(MWA) (ISalah et al.ll2005l ). The magnetic structure obtained from MWA measurements 
with this method will fill the missing link in coronal observations of the CME magnetic field 
and also place strong constraints on CME initiation theories. 



2. Modeling the Helios Observations 



The FR technique uses the fact that a linearly polarized radio wave propagating through 
a magnetized plasma will undergo a rotation in its plane of polarization. The rotation angle 
is given by Q = X 2 RM, where A is the wavelength of the radio wave. The rotation measure, 
RM, is expressed as 



RM 



37r 2 e m2c 3 



n P B • ds, 



where e is the electron charge, €q is the permittivity of free space, m e is the electron mass, c 
is the speed of light, n e is the electron density, B is the magnetic field, and <is is the vector 
incremental path defined to be positive toward the observer. FR responds to the mag- 
netic field, making it a useful tool to probe the coronal transient and quiet magnetic fields. 
Note that the polarization vector may undergo several rotations across the coronal plasma. 
Measurements at several frequencies are needed to break the degeneracy; observations as a 
function of time can also help to trace the rotation through its cycles. 

In situ observations of CMEs from interplanetary space i ndicate that CMEs are often 
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1990). This helical structure either exists bef ore the eruption (IChen 



19961 ; iGibson fc Lowlll998l ; iLin fc Forbesll2000l ). as needed for support 
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ing prominence mate rial, or is produced by magnetic reconnection during the eruption (e.g. 



Mikic &: Linker! 1 1 9941) . The flux rope co nfiguration reproduces the white-light appearance of 



CMEs f lChen 



1996; 



Gibson &: Lowlll998l ). This well-organized structure will display a specific 



FR signature easily discernible from the ambient medium, but direct proof of the flux-rope 
geometry of CMEs at the Sun has been lacking. 



2.1. Force- Free Flux Ropes 

Here we model th e Helios observations using a cylindrically symmetric force-free flux 



rope (ILundquistl Il950l ) with 

B = B J (ar)z + BoHJ^ar)} (2) 

in axis-centered cylindrical coordinates (f , 0, z) in terms of the zeroth and first order Bessel 
functions J an d J\ respectively, where B is the field magnitude at the rope axis, r is the 
radial distance from the axis, and H specifies the left-handed (-1) or right-handed (+1) 
helicity. We take ar = 2.405, the first root of the J function, so a determines the scale of 
the flux-rope radius r$. The electron density is obtained by assuming a plasma beta /5 = 0.1 
and temperatu r e T = 10 5 K, as implied by the extrapolation of in situ measurements (e.g., 



Liu et al.ll2005l . l2006al ). Combining equations (1) and (2) with a radio wave path gives the 



FR. 

For simplicity, we consider a frame with the x-y plane aligned with the flux-rope cross 
section at its center and the z axis along the axial field. Figure 2 shows the diagram of the 
flux rope with the projected line of sight. The flux rope, initially at 4i? away from the Sun 
with a constant radius r = 3.6i? and length 20i? Q , moves with a speed v = 500 km s -1 in 
the x direction across a radio ray path. The radio signal path makes an angle 9 with respect 
to the plane and <fi with the motion direction when projected onto the plane. The magnetic 
field strength at the rope axis is adopted to be B = 25 mG, well within the range estimated 



from the Helios observations (IBird et al.lll985l ) 



The resulting FR curves are displayed in Figure 3. A radio source occulted by the 
moving flux rope gives two basic types of FR curves, Gaussian-shaped and "N" -shaped (or 
inverted "N") depending on the orientation of the radio wave path with respect to the flux 
rope. When the radio signal path is roughly along the flux rope (say, for = 45° and 9 = 60° 
as shown in the right panel), the axial field overwhelms the azimuthal field along the signal 
path, so the FR curve would be Gaussian-like, indicative of a monopolar field. For a signal 
path generally perpendicular to the flux rope, the azimuthal field dominates and changes 
sign along the path, so the rotation curve would be "N" or inverted "N" shaped with a sign 
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change (left panel), suggestive of a dipolar field. These basic curves are consistent with the 
Helios measurements. Two adja cent flux ropes with evolving fields could yield a "W" -shaped 
curve as observed by Pioneer 6 (ILevy et al.lll969l ; iPatzold fc Birdlll998l ). The time scale and 
magnitude of the observed FR curves are also reproduced. When 9 = 0°, the line of sight is 
within the plane. Varying gives a variety of time scales of FR, ranging from ~ 3 to more 
than 10 hours, but the peak value of FR is fixed at ~ 57°. These numbers are consistent with 
the Helios data shown in the right panel of Figure 1. When 9 is close to 90°, the observer 
would be looking along the flux rope. The axial field produces a strong FR, but decreasing 
9 will diminish the rotation angle and make the curve more and more "N"-like. The time 
scale, however, remains at 4 hours. For <fi = 45° and 9 = 40°, the rotation angle is up to 
140°, in agreement with the Helios data shown in the left panel of Figure 1. 



2.2. Non- Force-Free Flux Ropes 



A non-force-free flux rope could give more flexibility in the field configuration. Consider 
a magnetic field that is uniform in the z direction in terms of rectangular coordinates. Since 
V • B = 0, the magnetic field can be expressed as 



B 



dA OA 

dy' dx' 



(3) 



where the yector potential is defined as A 
gives (e.g.. ISturrockill994l ) 



A(x, y)z. The MHD equilibrium, j x B — Vp = 0, 



d 2 A d 2 A 
dx 2 dy 2 



— A*o 



dA 



P 



2/i 



(4) 



where fio is the permeability of free space, p is the plasma thermal pressure, and j z is the z 
component of the current density. Equation (4) is known as the Grad-Shafranov equation. 
We see from this equation that p, B z and hence function of A alone. A special 

form of this equation V 2 A = exp(— 2 A) (in properly scaled units) has the solution (e.g., 



Schindler et al. 



1973|) 



A = \n 



a 



cos x + vi + a 2 cosh y 



(5) 



This nonlinear solution has been called the periodic pinch since it has the form of a 2-D 
neutral sheet perturbed by a periodic chain of magnetic islands centered in the current 
sheet. Here A, x and y are dimensionless quantities, and a is a free parameter that can be 
used to control the aspect ratio of the magnetic islands. 



- 6- 



From equations (3)- (5) we obtain 



B 



^oLq 



exp 



-2A 

BqLq 



B X = B C 



yl + a 2 sinh (y/L ) 



B n = B, 



acos(x/L ) + y/l + a 2 cosh(y/L ) 
a sm(x/L ) 



(6) 
(7) 



0" 



a cos(x I Lq) + y/l + a 2 cosh(y/L ) ' 
where Bq and Lq are scales of the field magnitude and length, respectively. The axial field 
B z and the thermal pressure can be obtained from 4r{p + ■§*-) = j z , which gives 



dA v' 1 2/^o ' 

2 r>2 / -2..n b\ 



B't B 2 
p + T~ = ^~ exp v r> / 



2/i 







where i?i is an arbitrary constant. Assuming a factor e in the partition of the total pressure, 
we have 



P 



Bl 



a cos ^- + Vl + ol 2 cosh y- ) + 



2/ 



B, 



-2 



a cos + Vl + ol 2 cosh -p- 



+ 



B[ 
Bl 



1/2 



(9) 



(10) 



Adjusting the parameters a and e gives a variety of flux rope configurations, circular and 
non-circular, force-free and non-force-free. 

A flux rope of this kind is displayed in Figure 4. As can be seen, this flux rope lies within 
a current sheet. To single out the flux rope, we require < x < 2ttLq and —ttLq/2 < y < 
7rLo/2 initially, where Lq = 1.5R Q . The flux rope is still 20R & long, moving with v = 500 
km s _1 across the line of sight. Other parameters are assumed to be Bq = 10 mG, B% = 0, 
a = 2, e = 0.1, and the temperature T = 10 5 K. Figure 5 shows the calculated FR. These 
curves are generally similar to those for a cylindrically symmetric force-free flux rope. Unlike 
the force-free flux-rope counterpart, the FR curves show a smooth transition from the zero 
angle to peak values. In addition, they are narrower in width, which may result from fields 
and densities which are more concentrated close to the axis. Note that the field magnitude 
is ~ 40 mG at the axis of the non-force-free flux rope. These profiles can also qualitatively 
explain the Helios observations. 

The above results suggest that CMEs at the Sun manifest as flux ropes, confirmin g 
what previously could only be inferred from in situ data (lBurlagalll988l ; lLepping et al.lll990l ). 
They also reinforce the connection of CMEs observed by coronagraphs with magnetic clouds 
identified from in situ measurements. 
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3. 2-D Mapping of CMEs 

As demonstrated above, even a single radio signal path can give hints on the magnetic 
structure of CMEs. Ambiguities in the flux rope orientation cannot be removed based on 
only one radio ray path. The power of the FR technique lies in having multiple radio sources, 
especially when a 2-D mapping of CMEs onto the sky is possible. 

3.1. A Single Flux Rope 

For a flux-rope configuration, the magnetic field is azimuthal close to the rope edge 
and purely axial at the axis. The rotation measure would be positive through the part of 
the rope with fields coming toward an observer and negative through the part with fields 
leaving the observer, so the azimuthal field orientation can be easily recognized with data 
from multiple lines of sight (radio ray paths). A key role is played by the axial component, 
which tells us the helicity of the flux rope. Consider a force-free flux rope for simplicity. For 
points on a line parallel to the rope axis within the flux rope, the field direction as well as 
the magnitude is the same. The fields on this line would make different angles with a variety 
of radio signal paths since the signal path is always toward the observer. As long as the 
axial field component is strong enough, these different angles will lead to a gradient in the 
rotation measure along the rope. 

Assuming an observer sitting at Earth, we calculate the FR pattern projected onto 
the sky for a force-free flux rope viewed from many radio sources. A flux rope has two 
possibilities for the axial field direction, with each one accompanied by either a left-handed 
or right-handed helicity. Plotted in Figure 6 are the four possible configurations as well as 
their rotation measure patterns. The angle 6 y defines the azimuthal angle of a line of sight 
with respect to the Sun-Earth (observer) direction in the solar ecliptic plane, while 9 Z is the 
elevation angle of the line of sight with respect to the ecliptic plane. The flux rope, with 
axis in the ecliptic plane and perpendicular to the Sun-Earth direction, is centered at 10i? Q 
from the Sun and has a radius of ro = 8R & and length 5Oi? . The magnetic field magnitude 
is assumed to be 10 mG at the rope axis. The gradient effect in the rotation measure along 
the flux rope is apparent in Figure 6 and it produces a one-to-one correspondence between 
the flux-rope configuration and the rotation measure pattern. The four configurations of a 
flux rope can thus be uniquely determined from the global behavior of the rotation measure, 
which gives the axial field orientation and the helicity. In order to fully resolve the flux 
rope, we have assumed ~ 80 radio sources per square degree on the sky, but in practice 
a resolution of 250 times lower can give enough information for the field orientation and 
helicity (see Figure 7). 
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The FR mapping obtained from multiple radio sources can also help to determine the 
speed and orientation of CMEs as they move away from the Sun. This mapping is similar 
to coronagraph observations. While the polarized brightness (Thomson-scattered, polarized 
component of the coronal brightness) is sensitive to the electron density, FR reacts to the 
magnetic field as well as the electron density and thus may be able to track CMEs to a larger 
distance than white light imaging. Figure 7 gives snapshots at different times of a tilted flux 
rope moving outward from the Sun. A Sun-centered coordinate system is defined such that 
the x axis extends from the Sun to Earth, the z axis is normal to and northward from the 
solar ecliptic plane, and the y axis lies in the ecliptic plane and completes the right handed 
set. A force-free flux rope, initially centered at (2, 2, 2)R Q in this frame and oriented at 30° 
from the ecliptic plane and 70° from the Sun-Earth line, moves at a speed 500 km s _1 from 
the Sun along the direction with elevation angle 10° and azimuthal angle 20°. The flux rope 
evolution is constructed by assuming a power law dependence with distance R (in units of 
AU) for the rope size and physical parameters, i.e., 

r = 0.2 x R - 78 AU 



B n = 15 x R~ 15 nT 



for the rope radius, 

for the field magnitude at the axis, and 

T = 3 x 10 4 x R-°- 72 K 

for the temperature. The rope length is kept at 3 times the rope diameter, and the plasma (3 
is kept at 0.1. Similar power-law dependences have been identified by a statistical study of 



CME evolution in the solar wind ([Liu et al.ll2005l l2006al ). but note that t he transverse siz e 



of the flux-rope cross section could be much larger than the radial width (ILiu et al.l l2006cf ) . 



The 2-D mapping has a pixel size of about 3.2 degrees. Even at such a low resolution, 
the flux rope can be recognized several hours after appearance at the Sun. The orientation 
of the flux rope with respect to the ecliptic plane is apparent in the first few snapshots, but 
note that this elevation angle may be falsified by the projection effect. The gradient effect 
in the rotation measure along the flux rope is discernable at 10 hours and becomes clearer 
around 20 hours. A right-handed helicity with axial fields skewed upward can be obtained 
from this gradient after a comparison with Figure 6 (top left). When the flux rope is closer 
to Earth, its appearance projected onto the sky becomes more and more deformed. Finally, 
when Earth is within the flux rope (around 80 hours), an observer would see two spots with 
opposite polarities produced by the ends of the flux rope. 



Note that the above conclusions are not restricted to cylindrically symmetric force-free 
flux ropes. We have also used the non-force-free solutions of the steady state Vlasov-Maxwell 
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equations (see §2.2), which unambiguously give the same picture. The FR technique takes 
advantage of an axial magnetic field coupled with the azimuthal component, which is the 
general geometry of a flux rope. This robust feature makes possible a precise determination 
of the CME field orientation. A curved flux rope with turbulent fields, however, may need 
caution in determining the axial field direction (see below). 



3.2. MHD Simulations with Background Heliosphere 



The above FR calculation does not take into account the background heliosphere. In this 
sense, the 2-D mapping may be considered as a difference imaging between the transient and 
background heliospheres. Here we use for the FR c alculation 3-D ideal MH D simulations of a 
CME propagating into a background heliosphere ([Manchester et al.l 12004 ). The simulations 
are performed using the so-called Block Adaptive Tree Solar- Wind Roe Upwind Scheme 
(BATS-R-US). A specific heating function is assumed to produce a global steady-state model 
of the corona that has high-latitude coronal holes (where fast winds come from) and a helmet 
streamer with a current sheet at the equator. A twisted flux rope with both ends anchored 
in the photosphere is then inserted into the helmet streamer. Removal of some plasma in 
the flux rope destabilizes the flux rope and launches a CME. The numerical simulation with 
adaptive mesh refinement captures the CME evolution from the solar corona to Earth. A 3-D 
view of the flux rope resulting from the simulations is displayed in Figure 8. The magnetic 
field, as represented by colored solid lines extending from the Sun, winds to form a helical 
structure within the simulated CME. The field has a strong toroidal (axial) component close 
to the axis but is nearly poloidal (azimuthal) at the surface of the rope. 

A fundamental problem in CME studies which remains to be resolved is whether CMEs 
are magnetically connected to the Sun as they propagate through interplanetary medium. 
Most theoretical modeling assumes a twisted flux rope with two ends anchored to the Sun 
(jChenl Il996l ; iKumar &: Rustl Il996l ; iGibson fc Low! Il998l ) . This scenario is suggested by ener- 
getic particles of solar origin observed within a magneti c cloud ( Kahler Sz Reamesl 1991 ). A n 



isolated plasmoid is also a possible structure for CMEs ( Vandas et al. 1993a ; Vandasl 1993b ). 



The FR mapping is capable of removing this ambiguity in that it can easily capture a flux- 
rope geometry bent toward the Sun. To show this capability, we calculate the FR mapping 
of the simulated CME in a background heliosphere. The MHD model gives a time series of 
data cubes of 3OOi? in length. We subtract the background from the rotation measure of 
the CME data to avoid possible effects brought about by the finite domain. Figure 9 shows 
the difference mapping of the rotation measure at a resolution of ~ 3.2 degrees when the 
CME propagates a day (~ 70R Q ) away from the Sun. The simulation data are rotated such 
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that the observer (projected onto the origin) can see the flux rope curved to the Sun. The 
coordinates, 9 y and 9 Z , are defined with respect to the observer. A flux rope extending back 
to the Sun is apparent in the difference image. The outer arc with positive rotation measures 
is formed by the azimuthal magnetic field pointing to the observer while the inner arc with 
negative rotation measures originates from the field with the opposite polarity. The rotation 
measure difference is positive near the Sun, which is due to a pre-existing negative rotation 
measure that becomes less negative after the CME eruption. 

A closer look at the image would also reveal asymmetric legs of the flux rope. This effect, 
indicative of a right-handed helicity, is created by the different view angles as described above. 
The nose of the flux rope does not show a clear gradient in the rotation measure because 
the view angles of this part are similar. In the case of the two legs directed to the observer, 
two spots with contrary magnetic polarities will be seen, so the curved geometry may also 
help to clarify the field helicity. 

4. Summary and Discussion 

We have presented a method to determine the magnetic field orientation of CMEs based 
on FR. Our FR calculations, either with a simple flux rope or global MHD modeling, demon- 
strate the exciting result that the CME field orientation can be obtained 2-3 days before 
CMEs arrive at Earth, substantially longer than the warning time achieved by local space- 
craft measurements at LI. 

The FR curves through the CME plasma observed by Helios can be reproduced by a 
flux rope moving across a radio signal path. Two basic FR profiles, Gaussian-shaped with 
a single polarity or "N"-like with polarity reversals, indicate the orientation of the flux rope 
with respect to the signal path. Force-free and non-force-free flux ropes generally give the 
same picture, except some trivial differences reflecting the field and density distributions 
within a flux rope. The FR calculation with a radio signal path, combined with the Helios 
observations, shows that CMEs at the Sun appear as flux ropes. 

2-D FR mapping of a flux rope using many radio sources gives the field orientation 
as well as the helicity. The orientation of azimuthal fields can be readily obtained since 
they yield rotation measures with opposite polarities. The axial component of the magnetic 
field creates a gradient in rotation measure along the flux rope, with which the flux rope 
configurations can be disentangled. Time-dependent FR mapping is also calculated for a 
tilted flux rope propagating away from the Sun. The orientation of the flux rope as a whole 
and its projected speed onto the sky can be determined from the snapshots of the flux rope 
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mapped in FR. We further compute the FR mapping for a curved flux rope moving into 
a background heliosphere obtained from 3-D ideal MHD simulations. It is shown that the 
FR mapping can resolve a CME curved back to the Sun in addition to the field orientation. 
Difference imaging is needed to remove the FR contribution from the background medium. 

The global FR map is a new technique for measuring the CME magnetic field. This 
method can determine the magnetic field orientation of CMEs without knowledge of the 
electron density. The electron density could be inferred from Thomson scattering measure- 
ments made by the SECCHI instrum ent (suite of w ide angle coronagraphs) on STEREO 



which has stereoscopic fields of view (lHowardll2000l ). With the joint measurements of the 



electron density, the magnetic field strength can be estimated. 

Note that the above results are a first-order attempt to predict what may be seen in 
FR. An actual CME likely shows a turbulent behavior and may have multiple structures 
along the line of sight; the rotation measure, an integral quantity along the line of sight, 
could display similar signatures for different structures. Therefore, interpretation of the FR 
measurements will be more complex than suggested here. However, having an instantaneous, 
global map of the rotation measure that evolves in time will be vastly superior to a time 
profile along a single line of sight, and comparison with coronagraph observations and actual 
measures of geoeffectiveness (e.g., the D st index) for a series of real events will eventually 
lead to the predictive capability proposed in this paper. 

The present results also pave the way for interpreting future F R observations of CMEs 



by large radio arr ays, particularly those operating at low frequencies (jOberoi fc Kasperll2004 



Salah et al.ll2005l ). The MWA - Low Frequency Demonstrator, specially designed for this 



purpose at 8 0-300 MHz, will fe ature wide fields of view, high sensitivity and multi-beaming 



capabilities (ISalah et al.l 120051 ). This array will be installed in Western Australia (26.4°S, 
117.3°E), a radio quiet region. It will spread out ~ 1.5 km in diameter, achieving ~ 8000 
m 2 of collecting area at 150 MHz and a field of view from 15° at 300 MHz to 50° at 80 
MHz. The point source sensitivity will be about 20 mJy for an integration time of 1 s. 
The array is expected to monitor ~ 300 background radio sources within 13° elongation 
(~ 5Oi? ) from the Sun, providing a sufficient spatial sampling of the inner heliosphere. In 
addition, this array will be able to capture a rotation measure of ~ 10~ 2 rad m -2 and thus 
is remarkably sensitive to the magnetic field. Science operations of the array will start in 
2009. Implementation of our method by such an array would imply a coming era when the 
impact of the solar storm on Earth can be predicted with small ambiguities. It could also fill 
the missing link in coronal observations of the CME magnetic field, thus providing strong 
constraints on CME initiation theories. 
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Fig. 1. — Time profiles of FR (bottom) and spectral broadening (top) of the Helios 2 signal 
during the CMEs of 23 October 1979 (left) and 24 October 1979 (right) recorded at the 
Madrid station DSS 63. The apparent solar offset of Helios 2 is given at the top. The dashed 
vertical line indicates the arrival time of the CME leading edge with uncertainties given by 
the width of the box "LE". Large deviations in F R following t he lea ding edge indicate the 
arrival of the CME's bright core. Reproduced from iBird et al.l (119851 ). 
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Fig. 2. — Schematic diagram of a force-free flux rope and the line of sight from a radio source 
to an observer projected onto the plane of the flux-rope cross section. The flux rope moves 
at a speed v across the line of sight which makes an angle <ft with the motion direction. 
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Fig. 3. — FR at A = 13 cm through the force-free flux rope as a function of time. Left is the 
rotation angle with 9 fixed to 0° and <fi = [10°, 20°, 30°, 90°], and right is the rotation angle 
with fixed to 45° and 9 = [0°, 20°, 40°, 60°]. 
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Fig. 4. — Same format as Figure 2, but for a non-force-free flux rope embedded in a current 
sheet. 
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Fig. 5. — Same format as Figure 3, but for crossings of the non-force-free flux rope. 
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Fig. 6. — Mapping of the rotation measure corresponding to the four configurations of a 
flux rope onto the sky. The color shading indicates the value of the rotation measure. The 
arrows show the directions of the azimuthal and axial magnetic fields, from which a left- 
handed (LH) or right-handed (RH) helicity is apparent. Each configuration of the flux rope 
has a distinct rotation measure pattern. 
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Fig. 7. — FR mapping of the whole sky at a resolution of ~ 3.2 degrees as a tilted flux rope 
moves away from the Sun. Note that the motion direction of the flux-rope center is not 
directly toward Earth. Values of the rotation measure for each panel are indicated by the 
color bar within the panel. Also shown is the time at the top for each snapshot. 
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Fig. 8. — A 3-D rendering of the CME magnetic field lines at 4.5 hours after initiation. The 
color shading indicates the field magnitude and the white sphere represents the Sun. 



-22 - 




Fig. 9. — Mapping of the rotation measure difference between the MHD simulation at 24 
hours and the steady state heliosphere. The two color bars indicate the logarithmic scale of 
the absolute value of the negative (-) and positive (+) rotation measure, respectively. 



